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Abstract 

Impact craters exist on solid surface planets, their satellites and 
many asteroids. The aim of this paper is to propose a theoretical ex- 
pression for the product pr^vf, where the three symbols denote the 
mass densityradius and speed of the impactor. The expression is de- 
rived using well known results of solid state physics, and it can be used 
in estimating parameters of impactors which have led to formation of 
craters on various solid bodies in the Solar System. 



1 Introduction 



Craters of various sizes have been observed on the terrestrial planets,their 
satellites and the major part of the asteroids. The study of these craters, and 
the resulting constraints on the related impactors has become a separate field 
of research in planetary science, see for example pQ. One of the fundamental 
questions concerning the impactors is what can be concluded about them by 
combining astronomical data with results of solid state physics. 

Recent theoretical work ( for example [2]) has shown that the application 
of elementary principles of solid state physics to this problem gives physically 
plausible results in reasonable agreement with those obtained by celestial 
mechanics. It was assumed in [2] that the material of the target was a 
crystal lattice, and the calculations were performed per unit volume. It was 
assumed there the condition for the formation of a crater is that the kinetic 
energy of a unit volume of the impactor has to be equal to or greater than the 
internal energy of a unit volume of the material of the target. The result was 
an expression for the minimal speed which an impactor of given parameters 
must have when hitting a target with a predefined set of parameters, in order 
to form a crater. 

The aim of the calculation reported here is to take into account the di- 
mensional effects - to consider both the impactor and the crater it forms as 
objects of finite dimensions. The novelty of the approach discussed in the 
next section, compared to existing work such as [3], is the generality - it is 
based on principles of solid state physics, and it can be applied to any solid 
material. Formation of impact craters is here discussed from the viewpoint 
of pure solid state physics as the following analogous problem: what kinetic 
energy of an impactor is needed to produce a hole of given dimensions in a 
material with a predefined set of parameters? 



2 Calculations 

The physical keypoint of this calculation is the idea that the kinetic energy 
of the impactor must be equal to or greater than the internal energy of some 
volume, V2, of the target. The kinetic energy of the impactor of mass m\ and 
speed V\ is 
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where 9 is the angle of the trajectory of the impactor immediately preceeding 
the impact with respect to the vertical at the point of impact. The internal 
energy of a volume V 2 of a solid is [4J 

E 2 = ±mty 2 (2) 
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all the symbols have their standard meanings and V is the speed of ellastic 
waves in the material of the target. One of the problems consists in know- 
ing the volume V2, which is, in first approximation,equal to the volume of the 
crater. It is known from theoretical work [5] and references given there that 
the form of a crater is determined by a combination of "gravity scaling" and 
"strength scaling" , where the term "strength" reffers to the material strength 
of the target. Explicite expressions for the form of craters exist in j3] and they 
are applicable to 4 material types. In order to simplify somewhat the calcu- 
lations, and get slightly more general results, it was assumed in the present 
work that craters have the shape of a half of a rotating elipsoid,with distinct 
semi axes, denoted by a, b and c. Physicallya and b denote the semi axes of 
the "opening" of the crater, while c is the depth. In that approximation,the 
volume of the crater is obviously given by 
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V 2 = -nabc (3) 
Inserting eq.(3) into eq.(2) leads to 

vr 3 {k B r 2 y 

E * = 15W (4) 

Inserting the definition of the mass density (pi = ^-) in eq.(l) and assuming 
further that the impactor has the form of a sphere of radius r 1; the following 
expression for the kinetic energy of the impactor is obtained: 

114 2 
Ei = -mxvl = - x -nrlpi(vcos9) 2 = --Kpir\{v cos 9) 2 (5) 

Equating expressions (5) and (4) and grouping terms, gives the following 
expression for the velocity which the impactor with mass density p\ and 
radius T\ must have in order to form a crater of depth c and semi axes of the 
"opening" a and b in a target with temperature T 2 

Pl r 1 (vcos9) = — ^ abc (6) 
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This expression can be reformulated as 



Equation (6) has the advantage of grouping known or measurable physical 
quantities on the right hand side, and those which are unknown on the left 
side. The application of eq.(7) demands the knowledge of the equation of 
state (EOS) of the material of the target. Generally speaking,regardless of 
its detailed form, an EOS can be expressed in the following analytical form 

p{p) = T.<-r (8) 

i=o Po 

where a; are some coefficients, p is the density at some pressure P and p is 
the density at pressure P. Limiting this developement to first order terms, 
it follows that 

S^s^ r 3 / 2 [1-3^-1 (9) 

dp po aipo 

A well known example of the EOS of a solid is the Birch- Murnaghan EOS 
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{l + (3/4)(^-4)[(p/p ) 2 / 3 -l]} (10) 

where B = — F(|^) T is the bulk modulus of the material and B' = (|^)t 
is its pressure derivative. The symbols p and p denote the mass density of 
a specimen under consideration at the initial value of the pressure Po and at 
some arbitrary value P. The final form of eq.(7) assuming that the material 
of the target can be described by the Birch Murnaghan EOS is 
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This expression would have a somewhat simpler form under the assumption 
that the bulk modulus is pressure independent, that is B' Q = 0. As B ^ 
the last expression does not tend to infinity for any real material. 

More refined results could be obtained by using analytic approximations 
to the Helmholtz free energy, from which all the thermodynamic potentials 
can be derived. This is,for example, the approach used in the ANEOS EOS 
[5J. However, ANEOS uses a file of up to 40 parameters for the characteri- 
zation of materials,which complicates its applications. 

3 An example 

As a test, the procedure disucssed in this letter was applied to the case of the 
Barringer crater in Arizona,a well known impact crater. It was assumed that 
V = 7km/ s, which is the mean value of the measured velocity of seismic waves 
in the Earth. Parts of the impactor have been found, and it is known to be 
an iron-nickel meteorite with p\ = 8000kg /m 3 . Taking that the temperature 
of the target is T 2 = 300K, inserting the known values of the diameter and 
depth of the crater (a = b = 1.2km, c = 0.18km) and assuming that the 
radius of the impactor is 50m, it follows that, for vertical impact, the impact 
velocity was around 13.7km/ s. These values are in good agreement with 
those existing in the literature jTJ. Repeating the same calculation for the 
impact angle of 9 = tc/4 leads to the value of the impact velocity of 19. 4km/ s. 
This value of the velocity is very close to the value obtained in [8], but their 
value of the radius of the impactor is r = 40m. 

4 Discussion 

In this paper we have presented a simple procedure which gives the possi- 
bility of estimating the value of the product of the density and radius of the 
impactor and the speed of impact: p\rf(v cos 8) 2 in terms of various material 
parameters of the target and the impactor. The expression for this product 
has been derived using basic principles of solid state physics, without any spe- 
cial assumption(s) about the materials. In (at least some of the) terrestrial 
applications the density of the impactor and the inclination of the trajectory 
can (in principle) be measured or estimated, which gives the possibility of 
estimating the impact speed if the value of T\ can somehow be estimated. 
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The shape of the craters has been approximated as a half of a rotating 
elipsoid,regardless of the composition of the material of the target. As a 
first approximation, the speed of the seismic waves has been "put in" in- 
stead of being calculated from the equation of state of the material of the 
target. Calculating this speed would have demanded the precise knowledge 
of the chemical composition of the target materials, and of B and B' for 
them. The impact of a projectile into a target leads to heating,and possi- 
bly melting and even vaporisation of the target around the impact point. 
The thermodynamic result of an impact depends on the heat capacity of 
the material of the target and on the kinetic energy of the impactor. The 
heat capacity can be measured, or theoretically estimated, assuming prior 
knowledge of the chemical composition of the target. Knowledge of the heat 
capacity is vital for estimates of the temperature to which the target heats in 
the impact. Therefore, the question is which impacts lead to melting and va- 
porisation of the material of the target, and which just provoke heating of the 
target. Even if the material of the target gets partially melted and as such 
flows away, the dimensions of the crater formed by the impact are determined 
by the volume of the material pushed away in the impact. This volume is, in 
turn, determined by the ratio of the kinetic energy of the impactor to the 
internal energy of some volume of the target. Some details of this problem 
have been discussed in [9]. 

The general conclusion of work reported here is that by using simple 
well known results of solid state physics and a given appoximation of the 
shape of impact craters, it becomes possible to estimate the value of the 
product p\r\ [v cos9) 2 . This in turn can be used to draw conclusions about 
the impactors which made various craters on the solid surfaces of objects in 
the planetary system. 
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6 Note 



This paper was submitted to Earth Moon and Planets. It was rejected by one 
of the editors, without being sent to a referee, on the grounds that the equa- 



tions are known and that the link http:/ /impact. ese.ic.ac.uk/ImpactEffects/ 
offers more data. 
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